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ABSTRACT 



Context. Chemically peculiar A type (Ap) stars are a subgroup of the CP2 stars which exhibit anomalous overabundances of numerous 
elements, e.g. Fe, Cr, Sr and rare earth elements. The pulsating subgroup of the Ap stars, the roAp stars, present ideal laboratories to 
observe and model pulsational signatures as well as the interplay of the pulsations with strong magnetic fields and vertical abundance 
gradients. 

Aims. Based on high resolution spectroscopic observations and observed stellar energy distributions we construct a self consistent 
model atmosphere, that accounts for modulations of the temperature-pressure structure caused by vertical abundance gradients, for 
the roAp star 10 Aquilae (HD 176232). We demonstrate that such an analysis can be used to determine precisely the fundamental 
atmospheric parameters required for pulsation modelling. 

Methods. Average abundances were derived for 56 species. For Mg, Si, Ca, Cr, Fe, Co, Sr, Pr, and Nd vertical stratification profiles 
were empirically derived using the DDAFit minimization routine together with the magnetic spectrum synthesis code Synthmag. 
Model atmospheres were computed with the LLmodels code which accounts for the individual abundances and stratification of 
chemical elements. 

Results. For the final model atmosphere Teff = 7550 K and log(g) = 3.8 were adopted. While Mg, Si, Co and Cr exhibit steep 
abundance gradients Ca, Fe and Sr showed much wider abundance gradients between logTsooo = -1.5 and 0.5. Elements Mg and Co 
were found to be the least stratified, while Ca and Sr showed strong depth variations in abundance of up to r; 6 dex. 
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1. Introduction 

. lOAql (HR7167, HD 176232, HIP93179) is one of the bright- 
' est rapidly oscillating Ap (roAp) stars. These stars exhibit high- 
overtone, low-degree, non-radial p-mode puls ations with peri- 
ods in the range of 6-24 minutes (Kurtz & Martinez L l2000t 
' lAlentiev et al.L |2012|) . The most distinquished characteristic of 
. roAp pulsations is the dependence of the radial velocity (RV) 
■ amplitude and phase on the particular chemical element/ion. 
After discovery of stratified abundance distribution in the at- 
mospheres of Ap stars it became possible to connect the ob- 
served pulsational charact eristics with the line depth formation 
in stratified atmospheres (i Rvabchikova et all, |2002|) . However, 
to get a self-consistent model atmosphere of an Ap star one 
needs to take into account anomalous abundances of most el- 
ements of the Mendelejev table as well as the stratified distri- 
bution of those elements that provide significant contribution to 
the continuum and line opacities. Recently, such modelling was 
performed for a few roAp stars: a Cir (Kochukhov et al., 2009), 
HD 24712 dShulvak et al.Ll2009h and HD 101065 - Przybylski' 
star dShulvak et al.Ll2010h . rSaio et all (l2010l) showed that the use 
of the self-consistent model atmopshere of HD 24712 allowed 
to provide better modelling of pulsations and to derive the theo- 
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retical distribution of RV amplitudes and phases over the stellar 
atmosphere, compatible with the observed RV pulsational char- 
acteristics. 

A first abundance analysis based on high resolution spectra 
was performed by Rvabchikova et al. (200Q). The mean mag- 
netic field modulus 1.2 kG was determined by comparison with 
magnetic spectrum synthesis calculations by iKochukhov et al.l 
(2002). Pulsations in lOAql were di scovered in a search fo r 
stars in the northern hemisphere by Heller & Kramer (1988|. 
Pulsat ions in spectral lines were detected by K ochukhov et aP 
d2002h . The star exhibits very sharp spectral lines which is typi- 
cal for very slow rotators. Being one of the brightest roAp stars, 
10 Aql was an excellent targ et for a phot ometric campaign with 
the space telescope MOST dHuber et all 12008). From their de- 
tailed analysis Huber et al. could derive a lower limit for the 
rotation period, f > 1 m onth, which was later confirmed by 
the study of RV pulsations dSachkov et all 1200 8V The authors 
performed detailed measurements of RV pulsations in the atmo- 
sphere of 10 Aql and found a phase jump in the RV pulsations 
derived from lines of the rare-earth ( REE) elements Nd and Pr. 
According to theoretical pred ictions (iKhomenko & Kochukhov! 
12009! Is ousa & CunhaLl20TTh the phase jumps may be associated 
with the so-called nodal regions. To get a proper pulsation model 
for any roAp star one needs to derive a depth-dependence of the 
pulsational characteristics. 
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The main goal of the present paper is to obtain the chemi- 
cal and atmospheric structure of lOAql and global parameters 
such as radius, mass and luminosity based on the available high 
quality spectroscopic, spectrophotometric and photometric ob- 
servations. The paper is organized as follows: spectroscopic ob- 
servations and data reduction are shortly presented in Sect.|2] the 
methods of the analysis are given in Sect. |3] The resulting abun- 
dances, element stratification and atmospheric structure calcu- 
lations are presented in Sect. |4] In Sect. |5] we discuss the im- 
portance of the detailed abundance and stratification analysis for 
atmospheric structure modelling. 

2. Observations 

As lOAql is a southern hemisphere object, the star was in- 
cluded in a larger survey for high resolution spectra of sharp- 
Uned Ap stars with the ESO-VLT UVES spectrograph, 68. D- 
0254(A) (PI Kochukhov). Dichroic standard settings DIC#2 cen- 
tered at 346+580 nm and 437+860 nm were used to obtain a 
wavelength coverage from 3030 to 10400 A with only a few 
minor gaps. The spectral resolution resulting from the chosen 
sUt width of 0.5 arcseconds is 7? 80 000 with a signal-to-noise 
ratio of 450 - 550 . Dat a reduction of UVES observations (e.g. 



iRvabchikova et all 12008 ) was performed with the ESO UVES 
pipeline ('Balleste r et al. , 1200 0). Additional continuum rectifica- 
tion was performed using a dedicated IDL programme. 



3. Methods 

3.1. Abundance Analysis 

As a first step in this investigation, a detailed abundance analy- 
sis was made for 10 Aql in order to derive input parameters for 
the first calculation of an individu a l mod el atmosphere with the 
LLmodels code by IShulvak et al.l (l2004t) . This code takes into 
account individual abundances for all elements and can addi- 
tionally accept vertical distributions of chemical elements as in- 
put. Abundances were derived from equivalent width measure- 
ments of 490 lines obtained with the IDL programme ROTATE, 
a graphical interface for comparison of observed and s ynthetic 
spect r a, and a m odified version of the WIDTH9 code tKuruczl 
Il993t iTsvmbali IT99d) . which includes additional magnetic 
line broadening. Abundances for 56 species were calculated. 
Atomic line data were extracted from the Vienna Atomi c Line 
Databa se (VALD II. jKupka et alJ jl999l) : IRvabchikova et all 
d 1 9991) : iPiskunov et all d 1 995h ) . For this initial abundance a naly - 
sis atmospheric parameters from IRvabchikova et al.l(l2000b were 
used, i.e. Tf^ff- 7650 K, \ogg- 4.0 and a mean magnetic field 
modulus of < B >=1.2 kG. The results were then used to com- 
pute an individual model atmosphere. Comparison of observed 
Ha line profiles with the synthetic profiles based on this new 
model revealed the necessity to reduce the effective temperature 
by 100 K in order to obtain a better fit. Table [1] lists the result- 
ing abundances of chemical elements. These abundances were 
then used to identify elements that appear to be vertically strati- 
fied in the atmosphere of the star Elements which revealed sys- 
tematic differences between abundances derived from lines of 
different ionization stages are expected to be inhomogeneously 
distributed in the stellar atmosphere. 



3.2. Stratification Analysis 

Vertical abundance stratification is considered to be responsi- 
ble for an impossibility to model cores and wings of strong 



lines with developed Stark wings with the same abundance, 
or as an impossibility to reproduce lines belonging to differ- 
ent ionization stages of the same element or high- and low- 
excitation fines wi t h a c hemically homogeneous atmosphere 
(IRvabchikova et all l2003h . All three effects were detected dur- 
ing the initial abundance analysis and nine elements were in- 
cluded in subsequent further stratification analysis. 

Starting from the T^ff - 7550 K, log(g) = 4.0 model at- 
mosphere with a peculiar but vertically homogeneous chemi- 
cal composition, stratification profiles of Mg, Si, Ca, Cr, Fe, 
Co, and Sr were co mputed with the DDAFit IDL-based code 
jKochukhovl l2007l) . This code iteratively fits a vertical abun- 
dance distribution to an unlimited set of observed spectral lines. 
The stratification profiles are defined by four parameters: upper 
atmospheric abundance, abundance in deep atmospheric layers, 
the vertical position of the abundance jump, and the width of the 
transition region where abundances change between the two val- 
ues. All four parameters are optimized simultaneously for one 
element at a time, but a fixed stratification of other elements can 
also be taken into account. 

Element Co was included in the model calculations only 
after a fi rst consistent stratified model had been derived by 
iNesvacill (12009) . as hyperfine structure data, which are nec- 
essary to synthesize observed Co line profiles, were not avail- 
able at the beginning of our study. Stratification of rare earth 
elements was derived using a semi-empirical fitting procedure 
based on observed equivalent widths under NLTE assumption. 
Therefore the shape of derived abundance profiles may devi- 
ate from a single-step fu ncti on. The method was des cribed in 
iMashonkina et alJ (l2009l) and lMashonkina et alJ (l2005l) . 

3.2.1 . Line Selection 

During this first stratification analysis large numbers of spectral 
lines were inspected most carefully in order to assure the em- 
ployment of the best suitable lines with the most accurate atomic 
parameters available. For the stratification analysis strong and 
weak lines of different ions are used to sample a wide range of 
atmospheric layers. Lines showing unusually broadened wings 
and relatively shallow cores, which were suspected to be a result 
of vertical element stratification, were also selected. 

A list of the atomic parameters used for further analysis 
is given for each element in Table [3] in the Online material. 
Atomic line parameters were extracted from VALD using the 
long format option for magnetic spectrum synthesis. The orig- 
inal sources of atomic data included in VALD are given in the 
corresponding table caption. 

Part of the line sets selected for stratific ation analysis in 
this w ork correspond to the ones pre sented by Kochukhov et al.l 
( 120061) and lRvabchikova et all ( l2008l) . The latter have studied si- 
multaneously the effects of peculiar Ca isotopic mixtures on line 
shapes and stratification profiles in a large number of Ap stars. 
In the present investigation the same Ca lines were used for the 
analysis, with the exception of lines which would be sensitive 
to the isotopic mixture. Therefore, ignoring this effect does not 
have a negative influence on the results of the present study. 

In order to estimate which atmospheric layers are sampled by 
the selected set of lines, formation depths of each line were cal - 
culated based on the method proposed bv iAchmad et al.l(| 19911) . 
which was implemented in our spectrum synthesis code. Based 
on the calculated contribution functions, a range of atmospheric 
depths which contribute most to the opacities at given wave- 
length points was derived for each line profile used in our strat- 
ification analysis. Figure [1] shows how these so called line for- 
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mation depths change when Fe or Sr stratification are introduced 
in model atmosphere calculations. In the case of Fe a few lines 
were available in the high energy region to sample the upper at- 
mosphere layers and constrain the upper abundance value. High 
excitation Fe lines are indicated to be useful for sampling lower 
atmospheric layers. For comparison, a similar plot is shown for a 
typical Sr stratification. The available observed Sr lines appear to 
be sensitive only to a narrow atmosphere depth region in the ho- 
mogeneous as well as in the stratified case. Therefore the derived 
location of the jump and the deeper atmospheric abundance are 
considered to be more reliable than the upper atmospheric abun- 
dance, which is determined by the depth of the line cores of the 
strongest lines. 

For iron, a large number of observed spectral lines was avail- 
able. Lines from UV, visual and near infrared wavelength regions 
were selected to ensure a better sampling of the atmospheric lay- 
ers where the abundance jump occurs, and put better constraints 
on the upper abundance value. In order to test the effect of the 
choice of line list on the resulting stratification profile two Fe 
distributions were calculated from a complete line list and one 
excluding UV and infrared lines. The difference between lower 
atmospheric abundances, and locations and slopes of the abun- 
dance jumps was found to be small, i.e. within the formal errors 
of the fitting routine. The upper atmospheric abundance derived 
from the line list including strong UV lines was found to be 
1 .7 dex lower than the one derived from the visual lines only. 
However, as only two lines in the sample are really sensitive 
to these layers, the uncertainty of the derived abundance is ex- 
pected to be of the order of 1 dex. 

In order to obtain the best possible sampling of atmosphere 
layers for all elements, the final stratification line lists aimed to 
include as many high excitation, UV and infrared lines as pos- 
sible, in addition to weak and strong lines from the visual range 
and lines that exhibit unusually broad wings and shallow cores 
due to vertical stratification. 

3.2.2. Influence of calculation setup on profile fitting 
procedure 

For several cases we find that the fitting routine DDAFit pro- 
duces results with large nominal errors for abundances or width 
of the abundance jump, especially when the atmospheric layers 
are probed by only a few spectral lines. Therefore small changes 
in the model atmosphere setup could have a large effect on the 
resulting parameters, which might be misinterpreted as being 
caused by real chemical gradients. A few tests were performed 
to investigate the behaviour of the code in more detail. 

In one experiment, a set of Fe stratification profiles was cal- 
culated based on the same model atmosphere but with different 
depth scales used to compute numerical derivatives in the mini- 
mization part of the code. Here we used optical depth log T5000 or 
column mass M. The resulting profiles showed no difference be- 
tween abundance values in deep atmosphere layers larger than 
O.I dex, and location as well as slope of the transition region 
were comparable. The upper atmospheric abundance differed by 
0.5 dex. This difference could be removed by using a model at- 
mosphere with 288 layers instead of 72 in DDAFit calculations 
in logTsooo. As the increased number of layers also increases 
overall calculation times by up to a few hours, for the rest of 
the analysis column mass was used as a depth scale. Finally, re- 
sults obtained this way are plotted on an optical depth scale. In 
the case of relatively wide jumps, this transformation results in 
rather smoothed shapes of the distribution profile instead of a 



simple two-step profile, due to the non-linear relation between 
T50()o and M. 

In another experiment, we computed different stratification 
profiles for seven elements using two model atmospheres with 
the same temperature, gravity and chemical composition, but us- 
ing 80 or 118 layers. The resulting profiles were qualitatively 
identical for most elements except for Ca, Fe and Sr In those 
three cases the model with higher depth resolution resulted in 
smaller abundances (-5, -1.5, and -2.5 dex) for upper atmo- 
spheric layers, while the overall distributions remained the same 
as for the model sampled on 80 depths. This finding indicates 
that for these three elements upper abundances are only con- 
strained by a few sensitive lines and small changes in atmosphere 
depth resolution might therefore have a large impact on the up- 
per atmospheric abundances derived with DDAFit. Adjusting 
the derivative step size used in DDAFit decreased the formal er- 
ror of the fitting procedure for the upper atmospheric abundance 
value below 0.8 dex. The difference of upper atmospheric Ca 
and Fe abundances between the two models thereby decreased 
as well from -1.4 to -0.7 dex, while it remained ^ -2.5 dex for 
Sr. This result indicates that the upper abundance of Sr is in- 
deed not very well constrained by the available set of spectral 
lines. Abundance values in higher atmospheric layers are very 
sensitive to small changes in model atmosphere and calculation 
setups for all three elements. Other fitting parameters for Ca, Fe 
and Sr, and overall stratification profiles of all other elements 
seem, however, not very much affected by such changes. 

During the iterative calculation of self consistent stratifica- 
tion models it is sometimes necessary to adjust T^^ to obtain bet- 
ter fits to observed fluxes. Similar to what was described above 
another test showed that a change of +100 K can result in a 
change of the upper atmospheric abundance by up to a few dex, 
if most available spectral lines are not very sensitive in these 
layers. Location and slope of the abundance gradient, as well as 
lower atmospheric abundances were not affected noticeably by 
such changes in T^g. 



3.2.3. Iterative Stratification Calculation 

The basic steps of the iterative procedure of atmospheric 
parameters and abu ndance determination were outlined in 
Shulvaketal. (2009). After the first set of stratification pro- 
files had been obtained, new model atmospheres for a set of 
[Teff , log(g)] pairs were calculated using the derived stratifica- 
tion as input. For testing the quality of various atmosphere mod- 
els observed hydrogen line profiles and spectral energy distribu- 
tions were used. In part icular we made use of spectrophoto metric 
observations by Adelm an et al.l (Il989l) and iBregerl (ll 9761) . ultra- 
violet energy distributions from the lUE satellite missiorfl, and 
spectrophotometric observations obtained by the STIS instru- 
ment mounted at the Hubble Space Telescope. A convergence 
of the stratification models was reached, when more iterations 
would not have improved the fit between synthetic line profiles 
and other observables any further, taking into account the formal 
error on abundances and step parameters of the DDAFit proce- 
dure. 



' http://archive.stsci.edu/iue/ 
" http://www.stsci.edu/hst/stis 
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Fig. 1. Range of atmosphere depths that contribute most to formation of selected line profiles, based on contribution functions 
according to Achmad et al. (1991) for Fe (left column) and Sr (right column). Upper panels show example stratification profile 
and formation depth ranges of selected lines in stratified model atmosphere, lower panels show formation depth ranges for lines in 
homogeneous model atmosphere. 



4. Results 

4.1. Abundance analysis 

Mean abundances were computed for a converged stratified 
model and are summarized in the last three columns of Table [1] 
For comparison, the results of the abundance analysis based on 
the initial homogeneous model (Teff = 7550 K and \og{g) = 4.0) 
and an intermediate stratified model with lower temperature are 
listed (columns 2 and 3 respectively). Including stratification in 
the model atmosphere calculation did not affect the overall abun- 
dances of non-stratified elements. The difference between ho- 
mogeneous and stratified models stayed within 0. 1 dex for most 
elements, which is well within the accuracy of the method of 
abundance determination by fitting equivalent widths. The com- 
parison of abundances from two stratified models with different 
temperatures illustrates that the mean abundances of homoge- 
neously distributed elements were not affected by the changes in 
model structure between different iteration steps of the stratifi- 
cation fitting procedure. 



The final converged model contains stratification profiles of 
nine elements, including REEs. Only for Pr and Nd a chemical 
gradient with increased abundances in upper atmospheric lay- 
ers was necessary to fit the observed spectral line profiles. All 
other elements indicate opposite gradients with concentration in 
deeper layers. A summary of the stratification profiles of all ele- 
ments in the final model is shown in Figure |2] Stratification pro- 
files did not change after convergence was reached, compared to 
the initial homogeneous model. 

The largest differences for upper atmospheric abundances, 
which were found to be the most variable of the four parameters 
during the iterative process, were observed for Ca and Fe. 

Detailed results of the DDAFit fitting routine, including for- 
mal errors for all four derived stratification parameters (upper or 
lower atmospheric abundance, position and width of the abun- 
dance jump in the column mass scale used for calculations) are 
presented in Table |2] 
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Fig. 2. Stratification profiles of the final model withFeff = 7550 
K and \ogg- 3.8 (solid lines) and initial homogeneous model 
(dotted lines). Dash-dotted lines indicate elements that were not 
included in the initial model but added at a later stage. 



Two examples of observed lines, synthetic profiles with the 
best fitting homogeneous abundance and stratified abundances 
are shown for Fe and Si in Fig. |3] Additional figures for other 
elements are included in the online material. Deviation between 
observed and final synthetic line profiles were between 1.7% and 
3.4% for all elements. 

4.1.1. Magnesium 

Magnesium stratification profiles did not change significantly 
between the first and final iterations and were not sensitive to 
small atmospheric structure changes in any region. Analysis of 
the contribution functions indicates that the line set is sensitive to 
a variety of atmospheric depths, as it includes two rather strong 
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Fig. 3. Fit of observed (black) and synthetic lines with homogeneous (red) and stratified (blue) abundances. Fe is shown in the left 
column, Si in the right column. 



Mg n lines in the UV as well as a number of strong infrared lines, 
in addition to strong lines in the optical spectral region. The final 
stratification profile shows an abundance jump of 1.7 dex which 
is needed to reproduce observed line profiles. 

4.1.2. Silicon 

A variety of suitable Si lines were available for stratification 
analysis including five high excitation Si ii lines useful for prob- 
ing deeper atmoshperic layers, two strong Si ii lines sensitive to 
upper atmospheric abundances, two infrared Si i lines as well 
as a number of Si i and Si ii lines that exhibit broad line wings. 
Throughout the iterative fitting process, the Si distribution re- 
mained stable between different iterations. The abundance jump 
of 2 dex is located in the same optical depth as for the Cr profile. 

4.1.3. Calcium 

Calcium shows the largest abundance stratification of all el- 
ements investigated in 10 Aql with an abundance change of 
5.7 dex between lower and upper atmospheric layers. The range 
of the transition region of the step profile was found to be 
logTsooo - [-0.2, -1.7]. The upper atmospheric abundance was 
very sensitive to small changes in atmospheric structure occur- 
ing during the iterative model calculations, before reaching a 
converged solution. As many strong Ca lines with anomalously 
strong wings and shallow line cores were included in the analy- 
sis, the lower boundary of the abundance jump was constrained 
with better accuracy than the upper atmospheric abundance. 

4.1.4. Chromium 

Like Fe and Si, a variety of spectral lines was available for 
Cr analysis, including lines with high excitation energies, and 
strong UV lines which were found to be very sensitive to high 
atmospheric regions up to log t - -2 according to the calculated 
contribution functions. The derived abundance profile remained 
very stable over all iterations and shows a small abundance gra- 
dient of 2 dex occuring in the same depth region as the Si abun- 
dance step. 

4.1.5. Iron 

Of all the investigated elements the highest number of suitable 
spectral lines for stratification calculations was found for Fe. 
Therefore the line list is expected to sample a large depth range 



of the stellar atmosphere. Two weak UV lines, two infrared lines, 
strong and weak lines of different ionization stages as well as 
high excitation Fe ii lines were included and a large final abun- 
dance gradient of 4.6 dex was derived. Inspection of the cal- 
culated contribution functions however revealed that only two 
strong lines were sensitive to higher atmospheric layers around 
log T - 2 while most other lines were mainly found to be formed 
between logr - -0.3 and -1.5 in a stratified atmosphere. The 
upper atmospheric abundance was therefore very sensitive to 
small changes in atmospheric structure between different iter- 
ation steps. Nevertheless the lower atmospheric abundance, the 
position, and slope of the abundance gradient did not change sig- 
nificantly between first and final iterations. 



4.1.6. Cobalt 

The stratification analysis of Co became possible with the avail- 
ability of better atomic line data, allowing to take into account 
hyperfine structure in spectrum synthesis. The coresponding data 
were taken from Pickering (19961) (Coi) and fBergemann et al] 
(2010) (Co ii). This is especially important as Co stratification is 
expected to be present (as indicated by the results of the abun- 
dance analysis) but not very strong. The final Co stratification 
was best represented by an abundance jump of 1.8 dex. This 
small abundance gradient (which was based on a short line list 
including strong UV lines and intermediate and weak optical 
lines) proved not to be very sensitive to small changes in model 
structure. 



4.1 .7. Strontium 

Formation depth calculations predicted that the available set of 
Sr lines would be sensitive to a very narrow atmospheric region. 
Upper atmospheric abundance values were therefore found to 
be highly dependent on small changes in model structure and 
choice of input parameters for the fitting routine, throughout the 
iterative stratification fitting procedure. The final stratification 
profile of the converged model, however, did not differ signifi- 
cantly from the initial one. A large abundance gradient of 4.7 dex 
was fit to the observed line profiles by our algorithm. However, 
as the upper atmospheric abundance is constrained by only a few 
points in the line profiles of the strongest lines, the real size of 
the abundance jump might be overestimated. 
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Table 1. Mean abundances and standard deviations derived from 
observed EQW with the starting model, an intermediate model 
with Teff = 7450 K and log(^) - 4.0, and the final model with 
Teff = 7550 K and \og(g) = 3.8. Column "#" denotes the number 
of lines measured for each ion. Stratified elements are marked by 
asterisks. 



Table 2. Stratification parameters for 7 elements based on the 
final model. Upper and lower atmospheric abundance values, 
abundance step position and width in column mass scale and 
deviation between observed and synthetic spectra are given. 
Formal errors of the DDAFit procedure are given. 
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4.1.8. Rare earth elements 

For Pr and Nd we performed a NLTE stra tification analysis as 
described in iMashonkina et al.l (l2005l 1200 9) using a trial-and- 
error method and the observed equivalent widths of the lines of 
the first and second ions, therefore the shape of the abundance 
profile may diff'er from a single-step function. Possible magnetic 
intensification for strong Nd m lines was approximated by pseu- 
domicroturbulent velocity of 1.0 km/s. The observed and calcu- 
lated equivalent widths of Nd lines are illustrated in Fig. |4] 

Pr and Nd are the only elements which are enriched in upper 
atmospheric layers between logr = -3.6 and -6.0 and there- 
fore required inclusion of inverted step profiles, compared to 
the other elements. The observed differences between upper and 



4.2. Fundamental and atmospheric parameters 

Our determination of atmospheric parameters is based on the 
model fit to the observed energy distributions calibrated to ab- 
solute units (see Sect. I3.2.3I I. We note a discrepancy between 
space- and ground-based observations as can be seen from 
Fig- El In particular, data taken from the STIS archive are sys- 
tematically lower than that of, e.g., broad-band spectrophotome- 
try of Adelman et al. ( 1989) and Breger (1976) as well as fluxes 
computed from Johnson and Geneva photometric systems. In 
the attempt to fit these datasets one additional assumption about 
atmospheric chemistry was made. As follows from the com- 
putations of particle diffusion, helium always sinks in subpho- 
tosheric layers of A-F stars (Michaud et al., 1979). The same 
He depletion was found later also in the optically thin layers 
dLeblanc & MoninLl2004 suggesting that the atmospheres of CP 
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Table 3. Fundamental parameters for He-normal and He-weak 
models with T^ff corresponding to the best fits to observed fluxes. 



He-normal 





7450 K 


7550 K 


R 


2.67 ± 0.06i?o 


2.59 + 0.06/?o 


M 


1.64 ± O.O8M0 


1.55 ±0.07 Mo 


L 


19.81 ±0.90Lo 


19.66 + O.9OL0 




He-weak 


Tcfi 


7400 K 


7450 K 


R 


2.72 ± 0.06Re 


2.67 + O.OGRa 


M 


1.70±0.08Mo 


1.64±0.08Mo 


L 


20.01 ± O.9OL0 


19.81 ±0.90Lo 



stars of spectral types A-F may well be He-deficient. Having no 
spectroscopic clues about the true He stratification, we thus con- 
sidered both He-normal (i.e. solar) and He-weak atmospheres 
with adopted He abundance oflog(He/H) - -4 (decreasing He 
abundance below this value has only marginal or no effect at all). 

In addition, using observed spectral energy distributions and 
stellar parallaxes allowed us to derive the radius of the star. The 
latter is obtained by minimizing the deviation between observed 
and predicted fluxes for a given Teff and log(^). The parallax 
n = 12.76 + 0.29 mas was taken from van Leeuwen (2007). 

Making use of STIS data resulted in T^g = 7450 K, log(g) = 
3.8 both for He-normal and He-weak models. Formally, the He- 
weak model with \og(g) = 3.6 provided a better fit with slightly 
lower compared to the model with \og{g) = 3.8, but the 
computed mass of the star is then found to be O.93M0 which 
is too small for an A-type star. Furthermore, the Paschen contin- 
uum, which is a good indicator of the atmospheric temperature, 
is clearly requiring a hotter temperature, which we find to be 
T^etf = 7550 K. The examples of the model predictions are shown 
on Fig. |5] The respective fundamental stellar parameters are: 
R = 2.53 + 0.06/^0, M = 1.47 + O.O7M0, L = 17.78 + O.85L0 for 
Teff = 7450 K model and/; = 2.46+0.06/^0, M = 1.39+O.O7M0, 
L = 17.73 + O.88L0 for the Tefr = 7550 K model respectively. 
Note that the error-bars results from the parallax accuracy only. 

As a next step we fitted the data obtained from ground-based 
observations. For this purpose we u sed spectrophotometric ob- 
servations by Adelman et al. (1989*) and extended them by NIR 
points (A > 7000A) from Breger ( 1976). We did so because ob- 
servations of Adelman cover the region of the Balmer jump and 
thus provide a homogeneous (in the sense of calibration) set of 
points in this important spectral region. In this particular case 
we find that He-normal models with Teff = 7450 - 7550 K, 
logig) = 3.8 and He- weak models with T^ff = 7400 - 7450 K, 
logig) = 3.8 provide the best fit. The comparison of the ob- 
served and predicted fluxes is shown in Fig. |6] The respective 
fundamental stellar parameters for the He-normal and He-weak 
settings are given in table [3] 

Figure Q illustrates a fit to the observed hydrogen Balmer 
line profiles for several models. It is seen that generally a model 
with Teff = 7450 K provides an optimal fit to all three Balmer 
lines. On the other hand H/3 and Hy are best fitted with a slightly 
cooler Tetf = 7400 K He-weak model. 

Taking into account a substantial scatter in the observed en- 
ergy distributions, the derived atmospheric parameters of the star 
are T,ff = 7500 + 50 K, log(g) = 3.8+0.1. There are two 
major spectral regions that play a crucial role in the determi- 
nation of atmospheric parameters: the Balmer jump (which con- 
trols best the value of surface gravity) and the Paschen contin- 
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Fig. 8. Comparison of stratification profiles of 7 elements com- 
puted with He-normal (solid lines) or He-weak (dashed lines) 
atmosphere models. 



uum (which controls best the effective temperature). Therefore, 
it is important to use observed fluxes from the same data source 
that cover these two regions. Different calibration schemes ap- 
plied by different observers could result in a systematic scatter 
in the derived model parameters. For instance, fluxes obtained 
by Breger lay slightly above the fluxes by Adelman in the re- 
gion of the Balmer jump (see, e.g., top right panel of Fig. |5). 
Using them entirely in model fitting together with lUE data or 
UV points of Adelman results in a hotter best-fit temperature of 
Teff - 7550 - 7600 K, compared to the one found by fitting data 
of Adelman extended with NIR points of Breger. Although this 
is not a large discrepancy, still one should be aware of similar 
systematics when analysing other stars for which no homoge- 
neous datasets are available. 

The difference between space- and ground-based observa- 
tions results in radius uncertainity of AR = Q.14Rq (assuming 
the same Teff for reference models). This is a rather small dis- 
crepancy, which is of the same order as th e errors provided by 
modem interferometry, as found by, e.g., iPerraut et al.l (UOllh 
for another bright Ap star y Equ. Still, with improving inter- 
ferometric accuracy, such a scatter between different observed 
spectrophotometric datasets could become important. 

To evaluate the influence of a change in model structure with 
modified He abundance on the derived stratification, two sets of 
profiles were computed using the He-weak and He-normal mod- 
els with Teff = 7550 K and \ogg = 3.8. Results of this compar- 
ison are shown in Fig.|8] A small systematic shift of abundance 
jumps towards lower atmospheric layers is observed when a He- 
weak composition is assumed. The quality of the fit of synthetic 
line profiles to observations was not found to be different be- 
tween the two models. 



5. Discussion and conclusions 

We have constructed a model atmosphere for 10 Aql using a self- 
consistent iterative approach to derive atmospheric abundances 
and vertical chemical gradients by fitting models to the observed 
spectral lines and wide-range energy distributions. 

By analyzing line profiles in a single spectrum we are able 
to derive somewhat mean vertical abundance gradients, i.e. gra- 
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Fig. 9. Temperature structure of models computed with homoge- 
neous and stratified abundance patterns. 



dients averaged over the visible surface of a star It is to note 
that modern Doppler Imaging techniques (that rely on the rota- 
tionally modulated variability of line profiles) provide us with 
two-dimensional images of stellar surfaces and reveal non- 
uniform horizontal distributions of chemical elements and mag- 
netic fields in atmospheres of CP stars. It is therefore natural to 
think of horizontal inhomogeneities as a result of variable verti- 
cal gradients, which depend on the local magnetic field geome- 
try, hydrodynamic flows, and possibly other processes (that we 
don't yet fully understand), the resulting force that pushes ions 
upwards or makes them diff'use downwards in the stellar atmo- 
sphere becomes variable over the stellar surface. That means that 
the positions of abundance jumps as well as their amplitudes are 
then the functions of local surface coordinates too. Only a self- 
consistent 3D mapping of stellar atmospheres has the potential 
of providing vertical and horizontal abundance gradients simul- 
taneously from the same observed data sets. Modem theoretical 
models of atomic diffusion indeed predict a rather strong sen- 
sitivity of element stratification profiles t o the geometry of the 
surface magnet i c field (see, for exarnple [Leblanc et al.L 
lAlecian & Still2010l:lstift & AlecianL l2()m 

However, the sharpness of spectral lines and lack of promi- 



nent line profile and magnetic fiel d variations imply a very 
long rotationa l perio d for 10 Aql (iRvabchikova et al.[ l2005t 
ISachkov et all l2008h . In this situation we are not able to use 
Doppler imaging to constrain surface distributions of chemi- 
cal eleme nts. Very weak long- term line profile variability de- 
tected by (ISachkov et al.L l2008h suggests that horizontal abun- 
dance gradients in this star are not large or that the star is viewed 
from nearly the same aspect angle during its rotation cycle. We 
have estimated that this weak variability would result in a mean 
abundance difference of <0.05 dex. Therefore, we do not expect 
stratification analysis to be significantly aff'ected by horizontal 
inhomogeneities. 

Our approach of using step-like assumption on derived abun- 
dance gradients has obvious limitations. First of all it is limited 
to a very narrow atmospheric region where a given set of atomic 
lines used in stratification analysis is formed. It naturally can- 
not predict abundances in optically very thin and thick layers 
which we simply do not see with our instruments. Secondly, 
it is impossible to restore any kind of complex vertical distri- 
bution whose shape may strongly deviate from a step-like as- 



sumption (see, for exa r nple Leblanc et al.L 120091: lAlecian & Stiftl 
l20I0HStift & AlecianLl20T2ir However, what is important for us 
is to see a general trend of element distributions as a function of 
depths, i.e. whether a given element is being brought upwards 
or downwards in the line forming region of a stellar atmosphere 
and how strong the resulting abundance gradient is. In spite of 
its limitations this approach is suitable of tracking any kind of 
systematic as applied to different stars with, e.g., diff'erent tem- 
peratures and surface magnetic fields. 

The detailed analysis of the high resolution spectra ranging 
from UV to infrared regions indicated that robustness of the in- 
ferred vertical stratification profiles is very much dependent on 
the sensitivity of the selected line set to a wide range of atmo- 
spheric depths. As seen, e.g. in our analysis of Cr stratification, 
inclusion of high- and low-excitation lines helps to constrain pa- 
rameters describing a chemical gradient best. However, upper 
atmospheric abundances are sometimes still only constrained 
by a few lines which are sensitive to high atmospheric layers. 
Therefore the value of the upper atmospheric abundance might 
vary between an initial homogeneous model atmosphere and a 
stratified model found after the iterative calculations have con- 
verged. 

The mean abundances for non-stratified elements changed 
by about 0.1 dex between the homogogeneous starting model 
and all other stratified models. Abundance variations between 
different stratified models were negligible. 

The most striking change between our homogeneous start- 
ing model and the final stratified model atmosphere is seen in a 
systematic shift of the abundance jumps towards higher atmo- 
spheric layers in the stratified case. 

Between the first and final models, derived stratification pro- 
files were used to calculate synthetic flux distributions to be 
compared with a large number of observations. This comparison 
was used for fine tuning of fundamental parameters of 10 Aql. 
We find that stratification itself has little effect on the synthetic 
energy distribution of the star compared to models computed 
with individual and homogeneous abundances. The overall fit 
and derived radii were almost identical. Therefore individual 
abundances appear to be the most important ingredient in the 
SED modelling, at least in the roAp temperature region. On the 
other hand, we do find that the observed hydrogen lines could 
not be reproduced with chemically homogeneous model atmo- 
spheres, as illustrated on Fig. Q where we also show two He- 
normal models with Teff = 7550 K, \og(g) = 3.8 computed with 
stratified and homogeneous abundances. Therefore, the inclu- 
sion of vertical element stratification has to be preferred when 
trying to establish an atmosphere model for a single Ap star in 
order to reproduce spectroscopic and SED observation simulta- 
neously. 

Similar to the case of HD 24712 (IShulvak et al.L l2009h we 
find a strong influence of REE opacity in the surface atmo- 
spheric structure of stratified models. In particular, strong abun- 
dance gradients of Pr and Nd shown in Fig. |2]lead to a heat- 
ing of plasma and appearance of the characteristic temperature 
jump. This is illustrated on Fig. |9] for He-normal final models. 
Interesting to note that such a temperature jump was first empir- 
ically predicted by Kochukhov et al. ( 2002) in order to fit narrow 
region of Balmer line profiles between wings and core (the so- 
called core-wing anomaly). Unfortunately, the temperature jump 
predicted by our models is located way too h igh in the atmo- 
sphere compared to what was suggested by Kochukhov et al.l 
(2002) (see their Fig. 3). We hope that a self consistent NLTE 
modelling of REE stratification will help to improve existing 
models and finally match the predicted position of the tempera- 
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lure jump. Last but not least, the incorporation of a n atmospheric 
model with an inverse temperature gradient from IShulvak et alJ 
(2009) allowed to i mprove pulsation models of HD 24712 as re- 
cently computed bv lSaio et alJ (l2010l) . All this suggests that the 
temperature structure of Ap stars is indeed very much different 
from the "canonical" T - t relation and that the derived temper- 
ature jumps are real. 

One of the important results of the present study is the de- 
crease of the log(^) value from 4.0 to 3.8 once the fitting of SED 
is introduced in the atmospheric analysis. This is caused by the 
incorporation of realistic chemistry in model atmosphere com- 
putations as well as the use of homogeneous spectral energy dis- 
tributions (i.e. STIS data in this particular case) which cover the 
whole Balmer jump region. A similar decrease of surface grav- 
ity was also requi red in t wo pr evious investigations of Ap stars 
a Cir (lKochukhov ~etan l2009h and HD 24712 (Shulvaket al., 
[2009). In the latter case, recent pulsation modelling by Sa io et alj 
t2010 ) also suggested that \og{g) of HD 24712 should be lower 
than what we usually find from photometric calibrations of Ap 
stars. From our analysis we conclude that a decrease of \og(g) 
is consistent with all observations if the appropriate abundance 
pattern is included in model calculations. On the other hand, de- 
creasing surface gravity does not dramatically change stratifica- 
tion profiles (see our first and final profiles in Fig. |2)- 

One of the open questions that still needs to be answered is 
the true helium content in the atmospheres of Ap stars. Using 
available observational material it is impossible to distinguish 
between He-normal and He-weak models with high accuracy. 
Surface temperatures of roAp stars are too low to see He lines. 
The helium concentration formally influences the quality of the 
fit to energy distributions, but once plotted against low and mod- 
erate resolution data the difference is difficult to see, except for 
the Balmer jump region. In addition, hydrogen line profiles are 
also only marginally affected by He depletion. Pulsation calcula- 
tions may therefore assume He depletion and still be compatible 
with all available observations (spectroscopy, spectrophotome- 
try). Again, improved self-consistent diffusion models are prob- 
ably the right way to go, in the search for a definite answer 

Acknowledgements. Part of this work was supported by the Austrian Fonds zur 
Forderung wissenschaftHcher Forschung, project P17890 and P22691. DS is 
granted by Deutsche Forschungsgeineinschaft (DFG) Research Grant RE 1664/7- 
1 . OK is a Royal Swedish Academy of Sciences Research Fellow, supported by 
the grants from Knut and Alice Wallenberg Foundation and Swedish Research 
Council. 

We also acknowledge the use of electronic databases (VALD, SIMBAD, 
NASA's ADS) and cluster facilities at Vienna Institute for Astronomy and Georg 
August University Gottingen. 



References 

Achmad, L., de Jager, C. & Nieuwenhuijzen, H. 1991, A&A250, 445 
Adelman, S. J., Pyper D. M., Shore, S. N., White, R. E., & Wanen, W. H., Jr 

1989. A&AS, 81,221 
Alecian, G., &; Stift, M. J. 2010, A&A, 516, A53 

Alentiev, D., Kochukhov, O., Ryabchikova, T., Cunha, M., Tsymbal, V., & and 

Weiss, W. 2012, MNRAS, 421, L82 
Artru, M.-C. Jamar C, Petrini, D., & Praderie, F 1981, AAP, 96, 380 
Ballester, P, Modigliani, A., Boitquin, O., Cristiani, S., Hanuschik, R.,Kaufer, 

A., & Wolf S. 2000, The ESO Messenger 101, 31 
Bergemann, M., Pickering, J. C, & Gehren, T. 2010, MNRAS, 401, 1334 
Bikmaev, I. F., Ryabchikova, T. A., Bruntt, H., Musaev, F. A., Mashonkina, L. I., 

Belyakova, E. V., Shimansky, V. V., Barklem, P. S., Galazutdinov, G. 2002, 

Astronomy and Astrophysics 389, 537 
Breger, M. 1976, ApJS, 32,7 

Buder K., Mendoza, C. & Zeippen, C.J. 1993, J. Phys. B, 26, 4409 
Fuhr J.R., Martin G.A., & Wiese W.L. 1988, J. Phys. Chem. Ref Data 17, Suppl 
4 

Garcia G. & Campos J., 1988, JQSRT 39, 477 



Den Hartog, E. A., Lawler J. E., Sneden, C, Cowan, J. J., 2003, ApJS, 148, 543 
Heller, C. H. & Kramer K. S. 1988, PASP 100, 583 

Huber, D., Saio, H., Graberbauer M., Weiss, W. W., Rowe, J. F, Hareter M., 
Kallinger, T., Reegen, P., Matthews, J. M., Kuschnig, R., Guenther D. B., 
Moffat, A. F J., Rucinski, S., Sasselov, D., & Walker, G. A. H. 2008, A&A 
483, 239 

Ivarsson S., Litzen U., & Wahlgren M. 2001, Phys. Scr 64, 455 

Jonsson, G., Kroll, S., Persson, A., & Svanberg, S. 1984, Phys. Rev., A30, 2429 

Khomenko, E. & Kochukhov, O. 2009, ApJ, 704, 1218 

Kochukhov, O., Bagnulo, S., & Barklem, P S. 2002, ApJ, 578, L75 

Kochukhov, O., Landstreet, J. D., Ryabchikova, T., Weiss, W. W., & Kupka, F. 

2002, MNRAS 337, 1 
Kochukhov, O., Tsymbal, V., Ryabchikova, T., Makaganyk, V., & Bagnulo, S. 

2006, A&A 460, 831 
Kochukhov, O. P. 2007, Proceedings of the Intemational Conference, held in the 

Special Astrophysical Observatory of the Russian AS, Eds: I. I. Romanyuk 

and D. O. Kudryavtsev. p. 109-118 
Kochukhov, O., Shulyak, D., & Ryabchikova, T. 2009, A&A, 499, 851 
Kupka F, Piskunov N.E., Ryabchikova T.A., Stempels H.C. & Weiss W.W. 1999, 

A&AS 138, 119-133 
Kurtz, D. W., & Martinez, R 2000, Baltic Astron., 9, 253 
Kuracz, R.L., 1993, CDROM 18 
Kuracz, R.L. 1994, CDROM 20-22 
Lambert, D. L. & Luck, R. E. 1978, MNRAS 183,79 
Leblanc, F. & Monin, D. 2004, The A-Star Puzzle, 224, 193 
Leblanc, F, Monin, D., Hui-Bon-Hoa, A., &; Hauschildt, P H. 2009, A&A, 495, 

937 

Martin, G.A., Fuhr, J.R., & Wiese, W.L. 1988, J. Phys. Chein. Ref Data, 17, 
Suppl. 3 

Mashonkina, L., Ryabchikova, T, & Ryabtsev, A. 2005, A&A, 441, 309 
Mashonkina, L., Ryabchikova, T., Ryabtsev, A., & Kildiyarova, R. 2009, A&A, 
495, 297 

Michaud, G., Martel, A., Montmerle, T, et al. 1979, ApJ, 234, 206 

Nesvacil, N., 2009, Diffusion in atmospheres of magnetic stars. Doctoral Thesis 

at the University of Vienna 
Perraut, K., Brandao, I., Mourard, D., et al. 2011, A&A, 526, A89 
Pickering, J. C. 2006, ApJS, 107, 811 

Piskunov N.E., Kupka R, Ryabchikova T.A., Weiss W.W. & Jeffery C.S. 1995, 

A&AS 112, 525-535 
Raassen A.J.J. & Uylings PH.M., 1998, A&A, 340, 300 

Ralchenko, Yu., Kramida, A.E., Reader, J., & NIST ASD Team. 

2008, NIST Atomic Spectra Database (version 3.1.5), available: 

http://physics.nist.gov/asd3 NIST, Gaithersburg, MD. 
Ryabchikova, T. A., Piskunov, N. E.. Stempels, H. C, Kupka, F, & Weiss, W. 

W. 1999, Physica Scripta, Volume T83. Issue 1, 162 
Ryabchikova, T. A., Savanov, I. S., Hatzes, A. P., Weiss, W. W, & Handler G. 

2000, A&A 357, 981 
Ryabchikova, T, Piskunov, N., Kochukhov, O., Tsymbal, V., Mittermayer P., & 

Weiss, W. W. 2002, A&A 384, 545 
Ryabchikova, T, Wade, G. A., & LeBlanc. F 2003, Proceedings of the lAU 

Symposium No. 210, Astronomical Society of the Pacific, 301 
Ryabchikova, T; Wade, G. A.; Aurire, M.; Bagnulo, S.et al. 2005, A&A, 429, 

L55-L58 

Ryabchikova, T, Ryabtsev. A.. Kochukhov, O., & Bagnulo, S. 2006, A&A, 456, 
329-338 

Ryabchikova, T, Kochukhov, O., & Bagnulo, S. 2008, A&A 480, 8 1 1 
Sachkov, M., Kochukhov. O., Ryabchikova, T, et al. 2008, MNRAS, 389, 903 
Saio, H., Ryabchikova, T., & Sachkov, M. 2010, MNRAS, 403, 1729 
Salih S., Lawler J.E., & Whaling W. 1985, Phys. Rev. A31, 744 (Co II) 
Seaton, M.J., Mihalas, D., & Pradhan, A.K. 1994, MNRAS. 266, 805 
Shulyak, D., Tsymbal, V., Ryabchikova, T., Sttz, Ch., & Weiss, W. W. 2004, 
A&A 428, 993 

Shulyak, D., Ryabchikova, T, Mashonkina, L.. & Kochukhov, O. 2009, A&A, 
499, 879 

Shulyak, D., Ryabchikova, T., Kildiyarova, R., & Kochukhov, O. 2010, A&A, 
520, A88 

Smith, G. 1988, J. Phys. B. At. Mol. Opt. Phys., 21, 2827-2834 

Sinith, G., & Raggett & D.St.J. 1981, J. Phys., B14, 4015 

Smith, G. & Gallagher A. 1966, Phys. Rev., 145, 26 

Smith, G. & O'Neil, J.A. 1975, A&A, 38, 1 

Sousa, J.C. & Cunha, M.S. 201 1, MNRAS, 414, 2576 

Stift, M. J., &; Alecian, G. 2012, MNRAS, 425, 2715 

Tsymbal V, 1996, in Astronomical Society of the Pacific Conference Series,Vol. 

108, M.A.S.S., Model Atmospheres and Spectrum Synthesis, S. J. Adelman, 

F Kupka, & W. W. Weiss, ed., pp. 198-199 
van Leeuwen, F 2007. A&A, 474, 653 

Wiese, W.L., Smith, W.M., & Miles, B.M. 1969, Atomic Transition Probabilities, 
vol. II, NSRDS-NBS 22, Washington, D.C. 



9 



N. Nesvacil et al.: A self consistent chemically stratified atmosphere model for the roAp star 10 Aquilae 
Wilke, R. 2003, PhD thesis, Heinrich-Heine-Universitat, Diisseldorf 



10 



N. Nesvacil et al.: A self consistent chemically stratified atmosphere model for the roAp star 10 Aquilae 
UV Balmer Jump 




2000 2500 
X, [A] 

Visuol and NIR 



3000 




3200 3400 3600 3800 4000 4200 
X, [A] 

IR 



t07450g3.80 

He-norm 

t07550g3.80 

He-norm 

t07450g3.80 

He-weok 

STIS 

lUE 

Breger 
TDl 

Adelmon 
Johnson 
Genevo 
2 MASS 



0.10 




5000 



6000 7000 8000 9000 10000 
X, [A] 



0.01 

000 1.01.21.41.61.82.02.2 
X, [/u.m] 



Fig. 5. Comparison between observed and predicted fluxes of 10 Aql. The model fit is done to the data obtained with STIS instru- 
ment. Please, see the online version of this figure for color-coding. 
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Fig. 6. Same as on Fig.|5]but with the model fit to the ground-based spectrophotometric data of lAdelman et al.l (Il989h and iBregerl 
([19761) . 
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Table 3. Atomic parameters for all lines used in stratification analysis of 10 Aql. Abbreviations of references correspond to different internal line 
lists of VALD and individual line lists with observed parameters. In case broadening constants were absent in VALD, approximations were used. 



Ion 


A [A] 


logg/ 


^low 


TStark 


Ref. 


Mgi 


4702.99 


-0.420 


4.346 


-4.460 


Bu 


Mgi 


4730.02 


-2.409 


4.346 


approx 


JK 


Mgi 


5172.68 


-0.380 


2.712 


-5.470 


Wi 


Mgi 


5183.60 


-0.160 


2.717 


-5.470 


Wi 


Mgi 


5528.40 


-0.400 


4.346 


-4.460 


Bu 


Mgi 


5711.08 


-1.833 


4.346 


approx 


Ku 


Mgi 


8806.75 


-0.200 


4.346 


approx 





Mgi 


8923.56 


-1.650 


5.394 


approx 


LL 


Mgii 


3104.71 


-0.030 


8.864 


-3.970 


Ku 


Mgii 


3104.72 


-1.330 


8.864 


-3.970 


Ku 


Mgii 


3104.80 


-0.190 


8.864 


-3.970 


Ku 


Mgii 


4390.51 


-1.700 


9.999 


-4.070 


Wi 


Mgii 


4390.57 


-0.530 


9.999 


-4.070 


Wi 


Mgii 


4427.99 


-1.200 


9.996 


-4.400 


Wi 


Mgii 


4433.98 


-0.900 


9.999 


-4.400 


Wi 


Mgii 


4481.12 


0.740 


8.864 


-4.700 


Ku 


Mgn 


4481.15 


-0.560 


8.864 


-4.700 


Ku 


Mgii 


4481.32 


0.590 


8.864 


-4.700 


Ku 


Mgii 


7877.05 


0.390 


9.996 


-4.540 


Wi 


Mgii 


8213.98 


-0.950 


9.999 


-4.770 





Si I 


5517.55 


-2.384 


5.080 


8.380 


BR 


Si I 


5645.03 


-2.400 


5.614 


approx 


NL 


Si I 


5665.56 


-2.040 


4.920 


8.290 


LL 


Si I 


5684.49 


-1.650 


4.950 


8.250 


NL 


Si I 


5690.43 


-1.870 


4.930 


8.300 


NL 


Si I 


5701.10 


-2.050 


4.930 


8.310 


NL 


Si I 


6125.02 


-1.513 


5.610 


approx 


sH 


Si I 


6131.57 


-1.705 


5.610 


approx 


sH 


Si I 


6131.85 


-1.685 


5.610 


approx 


sH 


Si I 


6142.48 


-1.420 


5.620 


approx 


BR 


Si I 


6155.13 


-0.770 


5.619 


approx 


BR 


Si I 


6244.47 


-1.363 


5.610 


approx 


sH 


Si I 


7017.66 


-1.070 


5.870 


approx 


sH 


Si I 


7034.90 


-0.780 


5.871 


8.140 


BR 


Sill 


5055.98 


0.440 


10.074 


9.040 


AJPP 


Sin 


5056.31 


-0.359 


10.074 


9.030 


NL 


Sin 


5632.96 


-0.820 


14.186 


approx 


Wil 


Sill 


5688.81 


0.000 


14.186 


approx 


Wil 


Sin 


5957.55 


-0.350 


10.067 


8.820 


AJPP 


Sin 


6347.10 


0.290 


8.121 


9.090 


Wil 


Sin 


6371.37 


-0.080 


8.121 


9.080 


NIST 


Sin 


6660.53 


0.230 


14.504 


approx 


Wil 


Cai 


4226.72 


0.244 


0.000 


-6.030 


SG 


Cai 


5857.45 


0.240 


2.933 


-5.420 


S 


Cai 


5867.56 


-1.570 


2.933 


-4.705 


s 


Cai 


6162.17 


-0.090 


1.899 


-5.320 


SN 


Cai 


6163.75 


-1.286 


2.521 


-4.998 


SR 


Cai 


6169.04 


-0.797 


2.523 


-4.997 


SR 


Cai 


6169.56 


-0.478 


2.526 


-4.994 


SR 


Cai 


6455.59 


-1.340 


2.523 


-6.072 


S 


Cai 


6471.66 


-0.686 


2.526 


-6.072 


SR 


K.cL II 




-i .zu / 


/ .J LJ 


/I Al 
-4.0iZ 


lr> 


Can 


6456.87 


0.410 


8.438 


-3.711 


TB 


Can 


8248.79 


0.556 


7.515 


-4.600 


TB 


Cri 


3355.73 


-1.107 


3.847 


-5.118 


GF 


Cri 


5065.91 


-1.340 


2.708 


-6.209 


GF 


Cri 


5110.74 


-1.289 


2.709 


-6.210 


GF 


Cri 


5204.50 


-0.208 


0.941 


-6.154 


MFW 


Cri 


5247.56 


-1.640 


0.961 


-6.120 


MFW 
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Table 3. continued. 



Ion 


D 

/I [A] 


log 5/ 


^low 


TStaik 


Ref. 


Cri 


5296.69 


-1.400 


0.983 


-6.120 


MFW 


Cri 


5345.80 


-0.980 


1.004 


-6.117 


MFW 


Cri 


5368.53 


-1.033 


3.847 


-6.104 


GF 


Cri 


6762.41 


-0.555 


5.280 


-5.246 


GF 


Cri 


6762.42 


-0.092 


5.282 


-5.259 


GF 


Cm 


3308.13 


-1.094 


4.920 


-6.627 


RU 


Cm 


3311.93 


-0.510 


4.156 


-6.648 


RU 


Cm 


3338.80 


0.053 


6.792 


-5.671 


RU 


Cm 


3355.65 


-0.809 


6.747 


-5.670 


RU 


Cm 


3378.33 


-1.045 


3.104 


-6.716 


RU 


Cm 


5024.52 


-1.733 


6.285 


-6.652 


RU 


Cm 


5085.65 


-2.215 


5.871 


-6.643 


RU 


Cm 


5210.82 


-3.259 


4.042 


-6.720 


RU 


Crn 


5210.86 


-2.941 


3.758 


-6.584 


RU 


Crn 


5232.49 


-2.360 


4.071 


-6.639 


RU 


Cm 


5379.79 


-2.081 


7.331 


-6.558 


RU 


Cm 


5678.39 


-1.496 


6.484 


-6.627 


RU 


Cm 


6418.90 


-1.915 


6.686 


-6.603 


RU 


Fei 


3354.05 


-1.169 


2.858 


-6.091 


V2 


Fei 


3356.68 


-1.554 


3.047 


-5.082 


GF 


Fei 


4404.75 


-0.142 


1.557 


-6.204 


V2 


Fei 


5198.71 


-2.135 


2.223 


-6.185 


V2 


Fei 


5324.17 


-0.103 


3.211 


-5.496 


V2 


Fe I 


5367.46 


0.443 


4.415 


-5.133 


V2 


Fe I 


5410.90 


0.398 


4.473 


-5.060 


V2 


Fe I 


5434.52 


-2.122 


1.011 


-6.303 


V2 


Fei 


5436.29 


-1.540 


4.386 


-4.997 


V2 


Fei 


5445.04 


-0.020 


4.386 


-4.582 


V2 


Fe I 


5546.50 


-1.310 


4.371 


-6.109 


V2 


Fe I 


5560.21 


-1.190 


4.434 


-4.323 


V2 


Fe I 


5576.08 


-1.000 


3.430 


-5.491 


V2 


Fei 


5862.35 


-0.058 


4.549 


-4.582 


GF 


Fei 


6136.61 


-1.400 


2.453 


-6.327 


V2 


Fe I 


6137.69 


-1.403 


2.588 


-6.112 


V2 


Fe I 


6165.36 


-1.474 


4.143 


-6.156 


V2 


Fe I 


6173.33 


-2.880 


2.223 


-6.194 


V2 


Fei 


6219.28 


-2.433 


2.198 


-6.202 


V2 


Fei 


6421.35 


-2.027 


2.279 


-6.310 


V2 


Fell 


4923.92 


-1.504 


2.891 


-6.583 


RU 


Fell 


5015.75 


-0.028 


10.348 


-5.287 


RU 


Fell 


5045.11 


-0.002 


10.308 


-4.984 


RU 


Fen 


5047.64 


-0.235 


10.308 


-4.976 


RU 


Fen 


5061.71 


0.284 


10.308 


-5.189 


RU 


Fen 


5197.57 


-2.348 


3.230 


-6.599 


RU 


Fen 


5291.66 


0.544 


10.480 


-5.468 


RU 


Fen 


5325.55 


-3.324 


3.221 


-6.603 


RU 


Fe n 


5362.86 


-2.616 


3.199 


-6.666 


RU 


Fe n 


5534.84 


-2.865 


3.245 


-6.601 


RU 


Fe n 


6432.68 


-3.687 


2.891 


-6.687 


RU 


Fen 


6516.08 


-3.432 


2.891 


-6.686 


RU 


Fen 


7449.33 


-3.488 


3.889 


-6.668 


RU 


Fen 


7711.72 


-2.683 


3.903 


-6.666 


RU 


Co I 


3502.27 


0.070 


0.432 


-6.300 


FMW 


Co I 


3502.61 


-1.240 


0.174 


-6.374 


FMW 


l_^0 1 








J.O^ / 


nF 

yjr 


Co I 


5342.69 


0.690 


4.021 


-4.928 


GF 


Co I 


5347.49 


-0.160 


4.149 


-4.952 


GF 


Co I 


6082.42 


-0.520 


3.514 


-5.474 


GF 


Con 


3415.77 


-1.740 


2.203 


-6.675 


SLW 


Con 


3501.71 


-0.970 


2.203 


-6.673 


T83av 


Sri 


4607.32 


-0.570 


0.000 


approx 


Bh 



14 



N. Nesvacil et al.: A self consistent chemically stratified atmosphere model for the roAp star 10 Aquilae 

Table 3. continued. 



Ion 


A [A] 




^low 


yStark 


Ref. 


Sri 


Mil. 11 


-0.220 


1.798 


approx 


GC 


Sri 


4741.91 


-0.410 


1.775 


approx 


GC 


Sri 


4784.31 


-0.510 


1.798 


approx 


Bh 


Sri 


4811.87 


0.190 


1.847 


approx 


GC 


Sri 


5111.19 


-0.380 


2.251 


approx 


GC 


Sri 


5229.26 


-0.330 


2.259 


approx 


GC 


Sri 


5450.83 


-0.340 


2.259 


approx 


GC 


Sri 


5486.13 


-0.460 


2.251 


approx 


GC 


Sri 


5504.17 


0.090 


2.259 


approx 


GC 


Sri 


5521.76 


-0.060 


2.251 


approx 


GC 


Sri 


5540.04 


-0.410 


2.259 


approx 


GC 


Sri 


6408.45 


0.510 


2.271 


approx 


GC 


Sri 


6503.99 


0.320 


2.259 


approx 


GC 


Sri 


lOlQ.Qil 


-0.030 


1.847 


approx 


GC 


Srii 


3380.70 


0.199 


2.940 


approx 


Bh 


Srii 


3464.45 


0.487 


3.040 


approx 


Bh 


Srii 


7334.95 


-1.102 


7.562 


approx 


Bh 


Ndii 


5130.59 


0.45 


1.304 


-5.61 


HLSC 


Nd II 


5165.13 


-0.74 


0.680 


approx 


HLSC 


Nd II 


5255.51 


-0.67 


0.205 


approx 


HLSC 


Ndii 


5293.16 


0.10 


0.823 


-5.74 


HLSC 


Ndii 


5319.82 


-0.14 


0.550 


-5.82 


HLSC 


Ndiii 


5102.16 


-0.62 


0.296 


approx 


RRKB 


Ndiii 


5294. 1 1 


-0.69 


0.000 


approx 


RRKB 


Ndiii 


5677.18 


-1.45 


0.631 


approx 


RRKB 


Ndiii 


5802.53 


-1.71 


0.296 


approx 


RRKB 


Ndm 


5845.02 


-1.18 


0.631 


approx 


RRKB 


Ndra 


5851.54 


-1.55 


0.460 


approx 


RRKB 


Ndm 


6145.07 


-1.33 


0.296 


approx 


RRKB 


Ndm 


6327.26 


-1.41 


0.141 


approx 


RRKB 


Ndm 


6550.23 


-1.49 


0.000 


approx 


RRKB 


Ndm 


6690.83 


-2.46 


0.460 


approx 


RRKB 


Prii 


4222.90 


0.271 


0.055 


approx 


ILW 


Prii 


5ill.ll 


-0.319 


0.483 


approx 


TT \1 T 

ILW 


Pr III 


5299.99 


-0.720 


0.359 


approx 


ISAN 


Priii 


6090.01 


-0.871 


0.359 


approx 


ISAN 


Prm 


6160.23 


-1.020 


0.173 


approx 


ISAN 


Prm 


6195.62 


-1.071 


0.000 


approx 


ISAN 


Prm 


7030.39 


-0.929 


0.359 


approx 


ISAN 


Prm 


7781.98 


-1.276 


0.000 


approx 


ISAN 







^ Stark data for Ca-Co are taken from Kumcz ( 1 9941) . Stark data for rare earth elements are tak en fro m iHartoe et al. I (l2003h . Abbreviations 
for references are as follows: Bu: ' Butler et al.i dl^sSTlK: iJonsson et alH ( fl984h ; Wi: |Wiese et al. ( 1969); Ku: 'Kumcz' (1993'); O: experimental 
comparison with solar spectrum and HD73666, T. Ryabchikova, L. Fossati (priv. comm.); LL: Lambert & Luck ( 1978); BR: Bikmaev et al. (2002,); 
NL: Kumcz (1993); sH: solar: Si=7.55 Holweger Model, E. Luck (private communication); AJPP: Artm et al. ( 1981); Wil: Wilke (2003); NIST: 
[Ralchenko et al. (2008); SG: Smith & Gallagher ( 1966); S: Smith ( 1988); SN: Smith & O'Neil ( 1975); TB: Seaton et al. ( 1994); SR: Smith et aQ 
(Il981i); G F: Kumcz ( 199 4); MFW: lMartin et al. ( 1988); RU: Raassen & Uvlings ( 1998); V2: various sources - see details in VALD publications 
iPiskunov et al.. (1995) ; Kupka et alJa999) ; .Rvabchikova et al. (1999); FMW : Fuhr et al. (.1988); SLW:ISalih et al I ( Il985h ; T83av: average valu e 
(Co 11) from CUNJ and SLW for the lin es with log(g/) > -1.0 ; Bh: TKumczl i\99t ): GC: lGarcia & Camposl ( 119881) ; HLSC: lHartog et al. 1 ( 120031) ; 
RRKB: [Ryabchikova et alJ t200d) ; lLW: llvarsson et all <200ll ); ISAN: Ryabtsev, private communication 
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Fig. 10. Fit of observed (black) and synthetic lines with homogeneous (red) and stratified (blue) abundances. Mg is shown in the left 
column, Ca in the right column. 





Fig. 11. Fit of observed (black) and synthetic Unes with homogeneous (red) and stratified (blue) abundances. Cr is shown in the left 
column, Co in the right column. 
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Fig. 12. Fit of observed (black) and synthetic Sr lines with homogeneous (red) and stratified (blue) abundances. 
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